Chiral Logarithms Tamed 
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We derive non- linear recursion relations for the leading chiral logarithms (LLs). These relations 
not only provide a very efficient method of computation of LLs (e.g. the 33-loop contribution is 
calculated in a dozen of seconds on a PC) but also equip us with a powerful tool for the summation 
of the LLs. Our method is not limited to the chiral perturbation theory only, it is pertinent for any 
non-renormalizable effective field theory such as, for instance, the theory of critical phenomena, the 
low-energy quantum gravity, etc. 

PACS numbers: 



The very fact of the spontaneous breakdown of the ap- 
proximate chiral symmetry in strong interactions leads 
to the possibility of the systematic expansion of hadronic 
amplitudes and correlation functions at low energies. The 
expansion is organized in powers of external momenta 
and the masses of pseudo-Goldstone bosons (pions). We 
denote the corresponding expansion parameter generi- 
cally as p 2 . An efficient method to perform the chiral 
expansion is based on the technique of effective chiral 
Lagrangian 1]. The leading 0(p 2 ) hadronic amplitudes 
can be obtained from the tree diagrams of the famous 
Weinberg Lagrangian Q: 



C 2 



tv[{d^Ud^) +m 2 {U + U^)] , 



(1) 



where F is the pion decay constant in the chiral limit, and 
m is the pion mass. The chiral corrections of the order 
0(p 4 In p 2 ) can be obtained from the one- loop calcula- 
tion with the Weinberg Lagrangian fT]) [l|, Q . In order to 
compute higher corrections, e.g. 0(p 4 ), 0(p e \np 2 ), etc., 
one has to include terms with four and higher deriva- 
tives in the effective chiral Lagrangian Note, how- 
ever, that the leading chiral logarithms (LLs), i.e. the 
correction of the form 0(p 2 [p 2 \np 2 ] n ) can be obtained 
from the n-loop diagrams generated by the Weinberg La- 
grangian The conciseness and beauty of the leading 
logarithm approximation lies in the fact that LL correc- 
tions depend only on one basic low-energy constant - 
F. In this approximation one avoids rapid proliferation 
of the low-energy constant with increasing of the chiral 
order. The calculation of LLs is a herculean task - it 
requires the computation of n-loop diagrams in the non- 
renormalizable field theory |T]). Presently, the LLs are 
computed to the two-loop accuracy for the 7T7r-scattering 
amplitude to the five-loop accuracy for the correla- 
tor of scalar currents Q , and to the three- loop accuracy 
for the generalized parton distributions (GPDs) [8|]. We 
note that for the case of the chiral corrections to GPDs 
the summation of LLs is indispensable 0, H| , because the 



smallness of the chiral expansion parameter is compen- 
sated by 1 /xgj . 

In this Letter we present an efficient method to com- 
pute LLs. This method reduces the task of LLs compu- 
tation to a simple algebraic problem, and paves a way 
towards the summation of LLs. 

We discuss details of the method for the massless 
0(N + 1) cr-model defined by the Lagrangian: 



C 2 = - [d^o-dao- + d^ a 8^ a ] , (2) 

where the fields are constrained by the relation a 2 + 
J2a=i ^ a ^ a = F 2 . We consider the 0{N + 1) cr-model 
for the following reasons: 

• It is equivalent to the massless two-flavour Wein- 
berg Lagrangian fTJ) for TV = 3. 

• It is a free field theory for N = 1, which can be 
used as a check of our calculations 

• It can be solved in the large- N limit, which provides 
a check of our calculations and allows us to asses 
the accuracy of the 1 /N expansion without tedious 
calculations. 

For simplicity, we consider LLs for the forward 7T7t scat- 
tering amplitude in the massless 0(N + 1) cr-model ([2]). 
The reader can easily apply our method to an observable 
in a field theory of her/his choice. The forward amplitude 
computed at Mandelstam t = has the form: 



T abca ( s ) = 5 ab 8 ca A{s) + S^S^Bis) + <J M <$ oc C(s). (3) 

The chiral expansion of the functions A(s), B(s), C(s) 
has the following structure 

oo n 

A(s) = (^) 2 S £5>«S"L fe , (4) 

oo n 

B(s) = (4tt) 2 S J2Y, B n )Sn Lk > C( S ) = A(-s), 



n=l k=0 
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where we introduce a dimensionless invariant energy 
S = s/(4ttF) 2 and L denotes the chiral logarithms, 
L = ln(/i 2 /s), /j is the renormalization scale. The 
first Aq -term in the expression for A(s) corresponds to 
the tree contribution to the scattering amplitude in the 
0(N + 1) a- model @. The other terms are higher chiral 
order corrections. Our aim consists in the calculation of 
the LL coefficients A„ = A n and B„ = B n appeared 
in Eq. (g|). 

Simple power counting [If shows that the coefficient 
An originates from the fc-loop diagram with vertices 
from the chiral Lagrangian C p with the number of deriva- 
tives of the pion fields p < 2(n + 1 — k). We see that the 
n-th order LL coefficient receives contribution from the 
n-loop diagrams with vertices generated by the leading 



Lagrangian ([2]). The UV divergencies in a n-loop dia- 
gram are removed by the subtraction of lower-loop graphs 
with insertion of the local counterterms corresponding to 
the subdivergencies of the original n-loop diagram. See 
detailed discussion of the structure of the subtractions 
m Refs. 0,0. The local counterterms relevant for our 
calculations renormalize the couplings the all-order La- 
grangian, which encodes the structure of counterterms: 



n=l 3=" 
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Here Pj are Legendre polynomials and a convenient no- 
tation for the operator monomials is introduced: 



J 



r,2n 



did 2 
d 2 



y = (tt^^ . . . aV Q ) d2(n ~ 3) ■ ■ ■ ^v b ) 



(6) 
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The coupling constants g n j (/j) are enumerated by two in- 
dices. The index n indicates the number of derivatives 
of the pion fields ( equal to 2n) in the corresponding 
countcrtcrm. We refer to the index n as "principal in- 
dex" . The second index j corresponds to the "exchanged 
spin" of the counterterm. The tree level contribution of 
the vertices ([3]) to the amplitude can be easily computed 
with the result: 



Atrcc( s ) 
■Btree(s) 



n=l o=o 
even 



(7) 
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even 
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The expansion coefficients A& and 15^' of the ampli- 
tude (|4]) are functions of the infinite set of couplings 
g n j(fj), denoted by g. These coefficients depend on the 
renormalization scale /i through /i-dependence of the cou- 
plings g. The renormalized (physical) amplitude is given 
by the sum of Eq. ^ and Eq. ([7]) and it must be in- 
dependent of /i. Thus, imposing the requirement that 
^ 2 ^( A ( S ) +A ct (s)) = we obtain (IH the following set 
of equations: 



n+l 



4 1) (g) + (-i) B E^yfe)=o 



3=0 
even 



(8) 



HA ( n k \ g ) + (k + l)A^ +1 \g)=Q. 

Here the /3-functions are defined as f3 n j(s) = 
H 2 -^p-g n j({i). Also we have introduced the differential 



operator H acting on the space of the coupling constants 

d 



n=l 3=o 
even 



(9) 



The set of equations (J8j> has the following solution: 



4 fc) (g) = 



(-1)' 



fe! 



The lowest constant A^ — gio — 1 is fixed by the tree 
level calculations with the Lagrangian ((2]). We see from 
the solution (TT0|) that, in order to obtain the LLs (con- 
stants A^ ) , we have to apply the operator H n times to 
a linear combination of the coupling constants g n j . This 
at first glance formidable problem can be solved if one 
notes the following crucial property of the operator H : 



-H> 



E 



9n+lj- 



(10) 



H n 9mj =0, if to < ft . 



(11) 



Indeed, the loop diagrams contributing to the renormal- 
ization of the constant g m j include vertices with con- 
stants g p i with p < to only. It implies that the (3 m j (g)- 
functions depend only on the subset of the low-energy 
constants g v \ with p < m. Hence, it is easy to see that 
the application of the operator H to g m j leads to the 
lowering of the principal index m by one unit, i.e. we 
can write the general form of the action of H 11 on g n +ij '■ 



H n g n+lj = n\ w 



nj9l0 ^' ^ n 



(12) 



The expression (|12[) simultaneously presents the defini- 
tion of quantities ui n j (n is the number of loops and 
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FIG. 1: One loop diagrams contributing to the /3-functions 
coefficients (|14p . Filled squares denote the counterterms (mi) 
and (pi) introduced in Eq. |(5j. 



functions, which is quadratic in couplings g n j, contribute 
to Eq. I|12p. This observation is in accordance with gen- 
eral RG analysis of Ref. Q . The general structure of the 
one-loop /3 n j-function is the following: 

n — 1 m n—Tn 

4 Tloop (g) = E E E B^ {n - m » 9mi9in - m)l (i3) 

m=l i=o 1=0 
even even 



j < (n + Mod(n, 2))), which determine the LL coefficients 
An = A { n } (see Eq. {ID]) with fc = n). 

Now our aim is to calculate ui n j from Eq. (fT2")l . 
Due to the property (fTTj) only one-loop piece of j3 n j- 



Here are numerical coefficients which can 

be obtained from the calculation of one-loop diagrams 
shown in Fig. 1 . The result of calculations gives the fol- 
lowing result: 



B 



2j+l 



N 



min[p,m] ifc ;fc fej 

_ (1 + 

fc=0 



2fc+ 1 



(14) 



where constants fi„ are computed as the following in- 
tegral with Legendre polynomials: 

fi n B = ^ T 1 j\dxP A (^jP B (x){ X -lY.(l$) 

The (n + 1) x (n + 1) matrix f2„ with matrix elements 
given by Eq. (I15|) . posses many beautiful and intriguing 
properties. For example, fi„ = 1 and tr(f2„) = sin(7rn), 
which follow from the fact that the matrix Q, represents 
SL(2, R) transformations. The relation of this symmetry 
group to the general structure of the renormalization pro- 
cedure in a wide class of effective quantum field theories 
will be discussed elsewhere. 

Now substituting the one-loop /3-function (fl~3f into 
Eq. (|12[) we obtain the following nonlinear recursive re- 
lations for desired coefficients Lo n j'- 

1 n — l m+1 n—m 

1 \ \ \ ~* „(m+l,i)(n—m,l) „s 
W ™J = ~ °) UmiU( n - m -x)lKW) 

n m=0 i=o 1=0 
even even 

The recursion (fT6|) allows to express the higher coef- 
ficients co n j through that with lower principal indices, 
starting with uiqq = 1. [We remind that n enumerates the 
loop order and j < (n+Mod(n, 2))] The coefficients of the 
/3-functions #( m +M)(n-m,z) are giyen by Eq Thc 

LLs for the amplitudes A(s),B(s) [(C n = (-l) n+1 A n )\ 
can be computed in terms of uj n j as follows: 

(2ti i o)| 



The recursive relation (|16p is the main result of the 
present paper. It allows a very fast computation of LLs. 
For example, the 33-loop chiral LL is computed in a 
dozen of seconds on a PC [l2|. The 6-loop results for 
LLs are presented [l3| in Table A for the amplitude A(s) 
and 7-loop results in Table B for the amplitude B(s). 

Our results for LLs for N = 3 agree with two-loop cal- 
culations of 7T7T amplitude || and with five-loop results 
for the correlator of the scalar currents [f|. Additional 
check of our method is provided by the case of N = 1. 
Indeed, for that value of N the Lagrangian ([2]) corre- 
sponds to a free field theory, therefore we should obtain 
nullification of LLs. For N—l the scattering amplitude 
is given by the sum A(s) + B(s) + C(s), it is easy to see 
that all LLs are cancelled in this case. 

The non-linear recursion relation (|16|) for LLs is valu- 
able not only because it renders a breakthrough in the 
calculations of LLs, it also allows general theoretical stud- 
ies of LLs anatomy in wide class of effective theories. As 
an example, let us consider the behaviour of LLs in the 
large- N limit of 0(N + 1) c-model. In this limit we can 
neglect all terms in the expression (fT4)) but the first one. 
Then, the recursion relation (| 16[) for large N LLs is sim- 
plified considerably: 

,LN 1 \ " JV , ,LN, LN /icf, 

~2^1^ 27TT a -W w («-™-iy (18j 

m— 1 

The solution of this recursion is obvious: 

<=(y)" V (19) 
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This solution is in agreement with the direct large- iV cal- 
culations in the 0(N + 1) cr-model, see e.g. In order 
to compute the 1 /N corrections to the leading result flT 
we substitute the LLs in the form: 



OjO 



c (1) 
N 



(20) 



into the recursion relation l|16p and obtain the linear 
equation for the coefficients c£). This equation has a 
solution in terms of the Lerch function. The correspond- 
ing expression is rather long, instead we give the result 
for the LL coefficients A n for the leading and subleading 
l/N orders in the case of large number of loops n ^> 1. 
The leading l/N asymptotic of the amplitude B we com- 
pute without any assumptions about n. The result is: 



N 

T 

'N 



, 7T 

¥ 



(i 



1 + 



It is a remarkable result! It shows that the l/N expansion 
for 0(N + 1) cr-model fails in the chiral order n ~ N and 
the expansion requires reordering. 



In summary, we have developed the method of non- 
linear recursion relations (|16[) which allows a calculation 
the leading chiral logarithms to essentially unlimited or- 
der. Furthermore, this method presents a puissant tool 
for study of general structure of infrared logarithms. It 
can be applied to any physics problem described by a 
non-renormalizable effective low-energy Lagrangian, e.g. 
theory of critical phenomena, low-energy quantum grav- 
ity, theory of magnetics, etc. 
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